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Abstract 

We investigate the profiles of diagonal components of static linear (efa and a w ), and first nonlinear (Jkxx and ffyf), optical 
responses of repulsive impurity doped quantum dots. The dopant impurity potential assumes Gaussian form. The dopant is 
considered to be propagating under damped condition which is otherwise linear inherently. The study principally puts 
emphasis on investigating the role of damping strength on the static polarizability components. In view of this the doped dot is 
exposed to a static electric field of given intensity. The damping strength interplays with the confinement potentials and 
fabricate the polarizability components delicately through the emergence of maximization, minimization, and saturation. Also, 
the difference in the extent of inherent confinement in x and i/-directions has been found to modulate the polarizability 
components. 
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Introduction 

Modem nanotechnology is advancing towards manufacture of new high-performance devices having significant 
technological applications. In view of this the linear and nonlinear optical properties of devices are found to be 
quite relevant. Low-dimensional quantum systems are well-known for exhibiting more prominent nonlinear 
optical effects than the bulk materials, and possess broad applications in high-speed electro-optical modulators, far 
infrared photo-detectors, semiconductor optical amplifiers and so on. Moreover, a lot of important information 
about the energy spectrum, the Fermi surface of electrons, and the value of electronic effective mass of these 
systems can be gathered on investigating their optical properties. In this context quantum dot (QD) has proved 
itself as a high performance semiconductor optoelectronic device. However, the properties of QD are strongly 
affected by the presence of impurities because of subtle interplay operating between QD confinement sources and 
impurity potentials. As a natural consequence we find a seemingly large number of important investigations in this 
field (Aichinger et al., 2006; Baskoutas et al., 2004; Giilveren et al., 2005; Movilla and Planelles, 2005; Rasanen et al., 
2004; Tas and Sahin, 2012). In practice, although the current technology is quite capable of manufacturing high 
purity QDs, the possibility of impurity contamination still persists. Particularly, during the fabrication of QD on 
the basis of colloidal solution, the impurity can infest the QD surface. 

The impurity potential modifies the dot confinement and thus engineers the electronic structure of these materials. 
In fact such modification can be achieved by tailoring the shape and size of the doped QD. The tailoring eventually 
changes the spatial disposition of energy levels and produces desirable optical transitions. Moreover, relevant to 
the optical transitions of QDs, the interplay between the dot confinement and impurity potential enhances the 
oscillator strength of electron-impurity excitation. Added to this, the optical transition energy depends on the 
confinement strength (i.e. the quantum size) and makes the resonance frequency more adjustable. These features 
are obvious for the manufacture of optoelectronic devices with controllable emission and transmission properties 
and ultra-narrow spectral linewidths. In consequence, impurity driven modulation of linear and nonlinear optical 
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properties assumes importance because of practical applications in photodetectors and high-speed electro-optical 
devices (Vahdani and Rezaei, 2009; Xie a , 2010). The fascinating optical properties of doped QDs have made them 
useful materials for various optoelectronic applications resulting in many important works on both linear and 
nonlinear optical properties of these structures (Qakir et al., 2010; Chen et al., 2013; Duque et al., 2013; Karabulut 
and Baskoutas, 2008; Khordad, 2010; Kumar et al., 2012; Sadeghi and Avazpour, 2011; Sahin, 2009; Vahdani and 
Rezaei, 2009; Yakar et al., 2010). 

The application of electric field provides pertinent information about the confined impurities (Lopez et al., 2008; 
Zeng et al., 2012). The electric field causes polarization of the carrier distribution with a consequent energy shift of 
the quantum states. The changed energy spectrum of the carrier could be harnessed to control and modulate the 
intensity output of optoelectronic devices. Furthermore, the electric field could destroy the symmetry of the system 
thereby leading to emergence of nonlinear optical properties. Among the nonlinear optical properties the second- 
order quantity is simplest and produces lowest-order nonlinear effect with magnitude usually stronger than 
higher-order ones, particularly if the quantum system demonstrates significant asymmetry. Thus, the applied 
electric field turns out to be extremely important so far as optical properties of doped QDs are concerned 
(Baskoutas et al., 2007; Duque et al., 2001; Karabulut and Baskoutas, 2008; Karabulut and Baskoutas, 2009; Lien and 
Trinh, 2001; Li and Xia, 2007; Murillo and Porras-Montenegro, 2000; Narayanan and Peter, 2012; Peter, 2006; Rezaei 
et al., 2011; Sadeghi, 2011; Vahdani and Rezaei, 2009; Xie and Xie, 2009; Xie a , 2010; Xie b , 2010; Zheng et al., 2012). 

The phenomenon of damping in QD bears a lot of importance in view of fundamental physics as well as 
nanoelectronic applications. The manufacture of high quality single electron transistors (Guo et al., 1997), logic 
elements (quantum bits) (Itakura and Tokura, 2003), memory cells (Yano et al., 1999), and lasers based on QD 
heterostructures (Dutta and Stroscio, 2000) is very much linked with damping. The phenomenon gains more 
importance since nanoelectronic devices are often introduced into large integral circuits consisting of densely 
packed structural elements (Fedorov et al., 2005). The typical length between the elements of integral circuits 
generally falls within the range of several tens of nanometers. These nanoelectronic devices mutually affect one 
another and are also influenced by the metallic and doped semiconductor fragments of the heterostructures and 
integrated circuits. It has also been observed that impurities are very much connected with damping (Li and 
Arakawa, 1997; Schroeter et al., 1996; Sersel, 1995) and different elementary excitation localized inside QD interplay 
with damping (Fedorov et al., 2005). The QD carriers have been found to be conspicuously interacting with 
environmental excitation in the vicinity of dopants. Consequently, we find a good number of studies on the QD 
electronic damping dynamics initiated by interference from environmental excitation (Itakura and Tokura, 
2003). There are also some experimental studies that examine damping in doped quantum confined structures by 
optical techniques (Ascazubi et al., 2002). 

In the present manuscript we have investigated the role of damping on the static linear and nonlinear 
polarizabilities of doped QD (Datta and Ghosh, 2011). The static optical response properties of these systems can be 
important in order to design devices with significant technological applications and can be principally useful in the 
area of all-optical signal processing. In the present work, for simplicity, we have considered inherently linear 
propagation of dopant from a fixed dot confinement center which is eventually impeded owing to damping. The 
damped propagation of dopant annihilates the inversion symmetry of the dot leading to the emergence of non-zero 
P value. The appropriate time-dependent modulation of spatial stretch of impurity (y _1 ) accompanying the 
damped movement has also been considered in this regard (Pal and Ghosh, 2013; Pal et al., 2013). The diagonal 
components of linear (a) and first nonlinear (p) polarizabilities have been computed by exposing the QD to an 
external static electric field. The present investigation principally addresses the role of damping strength (c) on the 
profiles of the direct components of linear ( a xx , a yy ) and first nonlinear ( P xxx , P yyy ) polarizability and reveals some 
important aspects. 

Method 

The model considers an electron subject to a harmonic confinement potential V (x, y ) and a perpendicular magnetic 
field B. The confinement potential reads V(x, y) = -m*o% (x 2 + y 2 ), where co 0 i s the harmonic confinement frequency. 
Under the effective mass approximation, the Hamiltonian of the system becomes 
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H ° = ^l~ ihV+ e -A] 2 +\m*a? 0 {x 2 + y 2 ). 


( 1 ) 


m* is the effective electronic mass within the lattice of the material considered. We have taken m* = 0.06 7m 0 and set 
tl=e = m 0 = a 0 = 1. This value of m*closely resembles GaAs quantum dots. In Landau gauge [A = (By, 0, 0)] ( A being 
the vector potential), the Hamiltonian reads 


H 0 = 


-±L(JL+ JL 


2m* \<5x 2 


dy‘ 


-) + ^m*of 0 x 2 + ^m*(of 0 + ®?)y 2 - t 


( 2 ) 


COc = — — being the cyclotron frequency (equivalent to magnetic confinement offered by B). The magnetic field in 

atomic unit corresponds to field strength of miliTesla (mT) order. We may define Q 2 = of + of as the effective 
frequency in the y-direction. The model Hamiltonian [cf. eqn (2)] represents a 2-d quantum dot with a single carrier 
electron (Chakraborty, 1999) with lateral confinement (parabolic) of the electrons in the x-y plane. The parabolic 
confinement potential often serves as an appropriate representative of the potential in semiconductor structures 
(Baskoutas et al., 2004; Giilveren et al., 2005; Lien and Trinh, 2001; Murillo and Porras-Montenegro, 2000; Peter, 
2006; Sadeghi and Avazpour, 2011; Yakar et al., 2010) and has been actually invoked in the study of optical 
properties of doped QDs (Qakir et al., 2010). The parabolic potential particularly becomes meaningful when the 
QDs are fabricated by etching process on a quantum well, by ion implantation or by application of electrostatic 
gates. 


In the present problem we have considered that the QD is doped with a repulsive Gaussian impurity (Halonen et 
al., 1996; Halonen et al., 1996). Incorporating the impurity potential to the Hamiltonian [cf. eqn (2)] we obtain 

H 0 ( x , y, ox, coo) =H 0 ( x , y, coc, coo) + Vim P (x o, yo) (3) 

Where Vimp (x 0 , yo) = Vimp (0) = Voe~ } 'of x ~ x o) + Cy _ yo) ] w jth y Q > 0 and Vo > 0 for repulsive impurity, and 
(x 0 ,y 0 ) denotes the coordinate of the impurity center. V 0 is a measure of the strength of impurity potential whereas 
y~f stands for the spatial stretch of the impurity potential. The presence of repulsive scatterer represents dopant 
with excess electrons. The use of such Gaussian impurity potential is quite frequent (Adamowski et al., 2005; 
Bednarek et al., 2003; Szafran et al., 2001). At this point it needs to be mentioned that Gharaati et al. (Gharaati and 
Khordad, 2010) introduced a new confinement potential for the spherical QD's called Modified Gaussian Potential 
(MGP) and revealed its effectiveness in determining the spectral energy and wave functions of a spherical quantum 
dot. 


Thus, the problem boils down to modeling the energy eigenvalues and eigenvectors of the two dimensional 
Hamiltonian H 0 : 

H 0 y/ n (x,y)= E n y/ n (x,y) (4) 

Equation (4) turns out to be the energy eigenvalues equation of a two dimensional harmonic oscillator as a> c (i.e. B) 
— > 0 and I q — > 0. Naturally, we seek diagonalization of // (l in the direct product basis of harmonic oscillator 
eigenfunctions. The time-independent Schrodinger equation has been solved using variational method expressing 
the trial wave function i/Xx, y) in the product basis of harmonic oscillator eigenfunctions (Datta and Ghosh, 2011; 
Pal and Ghosh, 2013; Pal et al., 2013) (f f ux) and <f> m (J3y) respectively, as 

yAx, y) = In.m Cn.m <!> n ( «*) <l m (Py) (5) 


Where C n m are the variational parameters and 


a 


lm r a> o 

I h 


and P = 


m*n 


. The general expression for the 


matrix elements of /// and V bnp (0) in the chosen basis are as follows (Datta and Ghosh, 2011; Pal and Ghosh, 2013; 
Pal et al., 2013): 


(^o)n,m; n m 

h j( n + ^)® 0 + (m'+ ^)Jo/ 0 + S nn S m m 
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And 


Where, 


And 


{V im J0)} , , = 

y J n,m;n ,m 

V 0 (<f> n {ax)<t, m m | e -r 0 [Cx-*o) 2 + (y-yo) 2 ] | ^ ax) ^ m ) 

= V 0 ( (j) n {ax)\e~ y o [C^-^o) 2 ] | </> n ,(ax)) 

= V oDl D 2 Yt n o {n ’ n hTo {m ’ m) f{k,n,n )g(l,m,m ) , 

( 7 ) 


/ ( \k,n,n ' ) : 


, n+n , 

2 k k\ n k C\C ( 1 - a* 2 )— k H n+n '- 2k (a lPl ) , 


S (U,m> 2*1! m iC m \c{l — /f 2 ) 


2 


- I 


H 


m+m' -21 


(P1P2) 


H n (x) Stands for the Hermite polynomials of ll 11 ' order. The various other terms read as follows: 


p = 


„ AA 1 n 2 BA 2 7i 2 

Al = — 7 — ,£>2 = — f , with A= 


B = 


1 j = exp 


( 2 m+m' m!m '! x y 

r 0 x o (4~ 4)1 


(2 n+n ' nln'lTr) 2 

- , = a 2 + x 0 , = 0- + x 0 . 


, 2 2 = exp 


y 0 yo (^2 jp) 


7 0 yo * a r? P 

= *’ (T= V % = 


( .2\V 2 ’ alld ( ,24/2 • 

(i-**) (l-/f ) 


In the linear variational calculation, a considerably large number of basis functions have been exploited after 
making the routine convergence test. 


The dopant propagation can be realized in terms of time-dependent dopant coordinates x 0 —> x 0 (t)and y 0 -> y 0 (t). 
Such a dopant propagation, in turn, makes y (the spatial spread of dopant) time-dependent (Pal and Ghosh, 2013; 
Pal et al., 2013) as the latter depends on the instantaneous dopant location and therefore reads as y(f) = 
y Q exp { — rjyjxQ (t) + y^ (t)J , where 77 is a very small parameter and y g is the initial value of y. The time- dependence 
thus indicates a lazy extension of the domain over which the influence of dopant is dispersed that ensues its drift. 
Now the time-dependent Hamiltonian reads 

H(t)=[H 0 - V tmp (0)]+ V^t), (8) 

Where 


p 1 (t) = p oe -Xt)[{*-*o(t)} 2 +{y-yo(t)} 2 ] (9) 

In the present work the intrinsic time-dependence of dopant propagation has been considered to be linear for 
convenient handling of the problem so that xo(t) = xo+at, yo(t) = yo+bt (Pal and Ghosh, 2013; Pal et al., 2013). The 
hindrance offered by damping to the dopant propagation has been modeled by introducing an exponentially 
decaying term phenomenologically through the parameter C, (a measure of damping strength) giving rise to x<>( t) = 
(xo+at) e^Pand yo(t) = (yo+bt)e~P . Such representation causes a faster fall of dopant propagation with increase in the 
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damping strength in comparison with the undamped motion. The spatial stretch of impurity (y _1 ) also changes 
accordingly. The matrix element involving any two arbitrary eigenstates p and q of Ho due to Vi(t) reads [for the 
expressions of various related terms see (Pal and Ghosh, 2013; Pal et al., 2013)]. 


v pT ( f ) = (v' p (x>y)\Vi(t)\r q (x,y)j 

~ Inm In'm' Q nm,p E n'm'.q 

{<! n (ax) | Vi(t) | (j) n , (ax) <j> m , (flyj) 

= V 0 (i)Y J )iX q (t) ( 10 ) 


The external static electric field Vo of strength e is now switched on with 

V 2 = - [s x \e\x + £ y \e\y] 

Where s x and s y are the field intensities along x and y directions. Now the time-dependent Hamiltonian reads 


The matrix elements due to Vi reads 


H(t) = [H 0 - V imp (0)] + V^t) + V 2 . 


(V2)n,m-.n',m' = ( V p (x, y) \v 2 \ yr q (x, y) 


( 11 ) 


( 12 ) 


r* c 

L nm,p o n'm',q 


II 

nm n'm' 

( (j) n ( ax) (p m (Py) \v 2 \ 0 n ,( ax) 0 , (fly) ) 


2nm 


r* c 

°i Z-tn'm' u nm,p L, n , m' ,q 


n'+l 


+ 



+ Sy Xnm 5] 


C* C 

n'm' L, n r m r ,q 


1 

~P 



& ; 


m'+l.m 


+ 



^n,n 


The time-dependent wave function can now be described by a superposition of the eigenstates of Ho , i.e. 

yAx.y.t) = I q a q (t) y/ q 


(13) 


and the time-dependent Schrodinger equation (TDSE) [cf. eqn (14)] containing the evolving wave function [cf. 
eqn(13)] has now been solved numerically by 6-th order Runge-Kutta- Fehlberg method with a time step size At = 
0.01 a.u. and numerical stability of the integrator has been verified. 


dw 

th — — — Hu/ 
dt v 


Or equivalently 


ihd q (t) = Ha q (t) 


(14) 


The numerical solution of TDSE [cf. eqn (14)] gives the time-dependent superposition coefficients with the initial 
conditions a P ( 0) = 1, a^O) = 0, for all q A p, where p may be the ground or any other excited states of Ho. The quantity 
P fc (t) = \a k (t)\ 2 can be used as a measure of population of k th state of Ho at time t. Actually P k (t) is the probability 
of observing the electron in the k th eigenstate and is therefore given by modulus square of the time-dependent 
superposition coefficients. 


The energy of the system E displays typical dynamics when represented as 

E(t) = I k E k (0)P k (t), 


(15) 


Where Ek(0) is the energy of k ,h eigenstate of Ho at t = 0. The time-average energy of the dot (E) in the presence of 
damping and static electric field is computed by numerically integrating the following expression: 
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<£■> = i/ o r (16) 

Where T is the final time up to which the dynamic evolution of the system is monitored. We have determined the 
values of (E(±s x ,s y = 0)), (E(±2s x , s y — 0)), {E(s x — 0 ,+Sy)), and {E(s x — 0, ±2£- y ))and used the data to compute 
the direct components of frequency dependent polarizability of the dot by the following relations obtained by 
numerical differentiation (Datta and Ghosh, 2011): 

a xx 4 = 5 ~(Em - \[{E{e x )> + <£(-£•*)>] + ±[(E( 2s x )) + {E{-2e x ))] , (17) 

And a similar expression for Otyytx^ . The diagonal components of first non-linear polarizability (the quadratic 
hyperpolarizability) are given by components that are calculated from the following expressions. 

Pxxx^x = 

[E(s x , 0) - E(-6 x , 0)] - \ [E(2s x , 0) - E(-2s x , 0)]. (18) 

And a similar expression is used for P yyy component. 

Results and Discussion 

Linear Polarizability (a): 

In this section we would like to discuss the influence of damping strength (g) on the direct components (a xx 
and a yy ) of linear polarizability. Fig. 1 displays the variation of a xx and a yy components with c . In the limit of 
extremely small damping strength ( g — > 0) both the components start from a very high value. However, both the 
components exhibit steep fall with further increase inp. At a damping strength of g ~ 10 -3 a.u. the linear 
polarizability components undergo saturation with further increase in damping strength. The initial sharp fall of a 
components with increase in g can be due to the fact that an increase in damping strength restricts the motion of 
dopant away from the dot confinement center. Consequently, the dopant suffers some enhanced confinement 
which reduces the dispersive nature of the system forcing the a components to diminish. Further increase in g 
impedes the dopant motion so severely that it is compelled to reside very close to the dot confinement origin. The 
strong dot-impurity repulsive interaction now increases the dispersive character of the system which may lead to 
an enhancement of a components. The saturation of a components beyond g ~ 10 _3 a.u. indicates kind of 
compromise between the reverse factors that govern the dispersive character of the system. 



FIG. 1: PLOT OF a COMPONENTS VS. g: (I) FOR a xx AND (II) FOR a yy . 

It needs to be further noted that the (/-direction suffers from inherently greater confinement in comparison with the 
x-direction [cf. eqn (2)]. As a result the a xx component displays higher magnitude than its [/-directional counterpart 
[%!]■ 

First Nonlinear Polarizability (fS): 

The following section describes the influence of damping strength on the direct components (P xxx and P yyy ) of first 
nonlinear polarizability. Figs. 2a and 2b reveal the pattern of variation of P xxx a n d P yyy as a function of g . Unlike 
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linear polarizability, here the components show some noticeable differences in their behaviors as g is varied. The 
inherently strong confinement in the y-direction appears to be the major cause for displaying such a behavior. The 
//^component evinces a minimization at a damping strength of g~ 2.0 x 10 4 a.u. [fig. 2a] whereas the P yyy 
maximizes at g~ 1.5 x 10 _4 a.u. [fig. 2b]. 
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FIG. 2: PLOT OF f) COMPONENTS VS. g\ (A) FOR p xxx AND (B) FOR p yyy . 

However, both the components ultimately settle to some steady value in the over damped domain. The observed 
behavior can be conceived as a resultant effect of damping strength, the confining forces and the dot-impurity 
repulsive interaction. At low damping strength the dopant remains considerably away from the confinement 
center because of its drift and feels little repulsive force. A large g value makes the drift highly hindered and the 
dopant remains close to the confinement center. The proximity of the dopant to the confinement center increases 
the symmetric nature of the system which could cause a huge fall in P xxx component leading to its minimization. 
On contrary, because of the said proximity the intense dot-impurity repulsive interaction now governs the P yyy 
component and maximizes it. The severity of the said repulsive interaction is also reflected by the enhanced 
magnitude of P yyy component (~ 10 5 a.u.) over that of P xxx (~ 10 4 a.u.). Thus, the emergence of maximization for 
Pyyy com P onen f indicates dominance of dot-impurity repulsive interaction over that of the confining forces. The 
saturation in the P components in the over damped regime manifests kind of balance between the several 
parameters that controls the optical responses. 

Conclusions 

The diagonal components of static linear and first nonlinear polarizabilities of impurity doped quantum dots have 
been investigated as the dopant propagates under damped condition. Special focus has been given on the 
dependence of the polarizability components on the damping strength. In practice, they are found to be dependent 
on the interplay between damping strength and the strength of the confining potentials. The damping strength 
basically arrests the motion of the dopant in the vicinity of the dot confinement center. This, however, depending 
upon the resultant effect of confining force and dot-impurity repulsive interaction, causes sharp fall, maximization, 
minimization, and saturation in the linear and nonlinear polarizability components. The intrinsically different 
extents of confinements in x and y-directions mingle with the restricted motion of the dopant because of damping 
and affect the polarizability. The results are thus quite interesting and expected to bear significance in the relevant 
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field of research. 
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